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CONSERVATION LAW FOR THE CAUCHY-NAVIER EQUATION OF 
ELASTODYNAMICS WAVE VIA FOURIER TRANSFORM 

NGUYEN VAN VINH AND NGUYEN TUAN MINH 



ABSTRACT. In this paper, we use the method of Fourier analysis to derive the for- 
mula of the total energy for the Cauchy problem of the Cauchy-Navier elastody- 
namics wave equation describing the motion of an isotropic elastic body The con- 
servation law of the total energy is obtained and consequently the global unique- 
ness of the solution to the problem is implied. 



Oh ■ 1. Introduction 

<. 

Methods of Fourier analysis are widely used for studying the theory of partial 

differential equations (PDEs). An important aspect is to find fundamental and clas- 
sical solutions of linear PDEs via Fourier transform. By the isometric property of 
Fourier transform, one can estimate energy functionals of a solution to some con- 
crete mathematical modelling problems in physics. Our main work is to derive the 
formula of the total energy for the Cauchy problem of the Cauchy-Navier elasto- 
dynamics wave equation describing the motion of an isotropic elastic body. For 
explanations about physical context related to the mathematical theory of linear 
CT) ', elasticity we refer the readers to I HI041 |Mus75[ ISom501 . 



Let us resume about the mathematical formulation and some necessary nota- 
tions. Remark that, the Newton's second law leads to the Cauchy's motion equa- 
tion of an elastic body, which takes the form 

V-cr + f = pu, (1.1) 

where cr = ((Tij(x, t)) ij= Y^ is the Cauchy stress tensor field, u = (ui(x, t)) i=T ^ is the 
displacement vector field, f = (fi(x, t)) i= j^ is the vector field of body force per unit 
volume and p is the mass density. 
The infinitesimal strain tensor field is given by the equation 

e = I [Vu + (Vu) T ] . (1.2) 

Moreover, the Hooke's law for homogeneous isotropic bodies has the form 

cr = A trace(e) I + 2/j, e, (1.3) 

where A,/i > are Lame's parameters and I is the second-order identity tensor. 
Substituting the strain-displacement equation (|1.2[) and the Hooke's equation (|1.3|) 
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into the equilibrium equation (|1.1|) , we obtain the Cauchy-Navier elastodynamic 
wave equation 

(A + ^)Vdiv(u) +^Au + f = pii. (1.4) 

This equation in the Cartesian coordinates has the form 



( A + / i )^(Etr)+^ + ^ 



0=1 



d 2 u k - — 



The Cauchy problem consists in finding a vector function u = [u\, u 2 , ..., u n ), which 
satisfies (|1.4|) and the initial condition as follows 



<9u 
u (M)lt= = ¥>o0), ■^■( x ' t ) 



¥>iW 



(1.5) 



i=0 



It is easy to show that the Cauchy problem (|1.4til.5|> has a Cauchy-Kovalevski- 
Somigliana solution in the following form 

d 2 



u 



X + 2p 



A w 



A + p 



Vdiv(w), 



dt 2 p j p 

where w is a solution to the Cauchy problem of the biwave equation 

(P_ _ A + 2p 
dt 2 ' 



(1.6) 



-A 



w 



t=0 



p 

dw 
~dt 



^-^A)w = fM 



0; 



(1.7) 
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w 



i=0 
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V»i- 



t=o 



Note that, in IIKVM111 , a formula of exact solution to the Cauchy problem of the 
biwave equation is given by using Fourier transform. 



The total energy E(t) of a solution to the equation (|1.4|) is defined as the sum- 
mation of the kinetic energy and the potential energy, where the kinetic energy 
functional is given by 



U(t) 



P 



d_ 

at' 



-u 



dx 



P 



£ 

fc=i 



du k 
dt 



dx. 



and the potential energy functional or the strain energy has the following form 

K(t) = / I -A trace (e) 2 + p trace (e 2 ) ) dx. 

We denote by <S> (M n ) the space of Schwartz functions. Remark that an indefinitely 
differentiable function <p is called Schwartz function when <p and all its derivatives 
are required to be rapidly decreasing, i.e. 



:,/3 = SUP 



.ret 



" ( iL)'« x) 
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for every multi-index a and (3. The Fourier transform of e iS(IR n ) is defined by 

?[</>] (£) = ?(£) = [ e~ l ^<j ) (x)dx. 



It is well-known by the Plancherel theorem that the Fourier transform is an isome- 
try in <S(R") with respect to L 2 -norm, i.e. 



(x) I dx 



10(0 1 2 ^, 



for every <\> G 5(M n ). The basic concepts of Fourier analysis and application for 
solving PDEs and study the conservation law of energy can be found in HFol95 
ISS031 



By using some Fourier analysis techniques, we will show the exact formula of 
the total energy, which depends only on the initial data and is independent of time. 
Consequently, the global uniqueness of the solution to the Cauchy problem of the 
Cauchy-Navier elastodynamics wave equation is proved. 

2. Main results 

Let us denote a 2 = (A + 2/i)/p, b 2 = fi/p and assume that f = 0, then (|1.4[) can be 
rewritten as the following homogeneous equation 

d 2 



dt 2 



- b 2 A u - (a 2 - 6 2 )Vdiv(u) = 0. 



(2.1) 



Note that, if u(x, t) e C 4 (M n x [0, oo)) is a solution to the equations (|2.1H1.5[) . then 
it is also a solution to the homogeneous biwave equation 



dt 2 



together the initial condition 



u 



*=0 



^0 



a 2 A 



du 



' dt 



dt 2 



v 



- b 2 A u = 



(2.2) 



i) 



i=0 



d 2 u 



dt 2 

d 3 u 



I dt 3 



i=0 



t=0 



cp 2 = (a 2 - b 2 ) V(div Vo ) + 6 2 VVo> 
(p 3 = (a 2 - b 2 ) V(div ip x ) + b 2 V 2 <p x . 



(2.3) 



Indeed, taking divergent to the both side of (|2.1[) , we have 

— div(u) - a 2 Adiv(u) = 0. 

Taking Laplacian and doubly differentiating with respect to t, we obtain that 

d 2 



dt 2 



:A - b l A z u 



^)VAdiv(u), 



(2.4) 



(2.5) 
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^-^)u=(^-^div(u). (2.6) 



From ((2ii l) - ((2i5l) - ((Z6)) , we have 



(&a - a 2 A 2 ] f JU - b 2 A 2 ) u 



V^ 2 / V^ 2 

= (a 2 - 6 2 )V fj^div(u) - a 2 Adiv(u)') = 0. 

The initial conditions (|2.3|) is easily verified from (|2.1H1.5|) . For more important 
discussions about the biwave equation, we refer to HKC07[ IKKC101 IKCT10H and 
IIKVM111 . 

Using the above connection, we will derive the formula of the total energy of a 
solution of (|2.1H1.5|) via Fourier transform. In the next sequence, we assume that 
u = (ui,...,Un), cp = (<po,i,...,<po, n ), Vi = (<Pi,i, ■-, <Pi,n) are vectors of Schwartz 
functions. 



Lemma 2.1. The Fourier transform of solution to the homogeneous biwave equation \2.2\ 
I2.3P has the following form 

~ ( t ,x COs(a|e|t) /fi r s, t , sin ( a KI^) (c ~ v t , 

u & «) = -^p— <e, vo) e + -^r- & vi> £+ 

i cos(6|£|t) , , , sin(6|£|t) , , v 

H |7p IKI Vo - (?> <A>) ?J H HTp — V&\ <Pi ~ (?> Vi) Sj • 

Proof. Taking Fourier transform to the both sides of the equation (|2.2|), we obtain 
that 

^u(e,t) + (a 2 + b 2 m 2 ^u(£,t) + a 2 b 2 \tfu(Z,t) = 0. 
This fourth ODE has the general solution, which takes the form 

u (f , t) = Ci cos (a\£\t) + C 2 sin (a\£\t) + C 3 cos (b\£\t) + C 4 sin (6|f |t) , 

where vector of parameters C\, C 2 , C 3 , C 4 are determined from the initial condi- 
tion (|2.3|) by solving a system of linear equations. After some simplifications, the 
Fourier transform of u is calculated as 

(2.8) 
In the above formula, we note that 

4>2 = -{a 2 -b 2 )(i,Cp )i-b 2 \i\ 2 ^ 
Substituting (|Z9]) into l[Z8]) , we get the formula (|Z7)) . D 
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Theorem 2.2. Let u be a solution o/ (12.2H1.5D . Then the total energy ofu takes the form 



E{t) =E(0) = I [ p\a 2 



i=l / i<j x J ' i=l / 



(2.10) 



Proof. We only need to verify that the total energy of a solution of the biwave equa- 
tion (|2.2H2.3|) also has the form (|2.10[) . Applying Plancherel theorem, the strain en- 



ergy functional can be rewritten as 

K(t) = / ( -A trace {if + /Urace(£ 2 ) j d£ 



We note that 



e 

2 



\ (£u T + uf) , e 2 = \ ((£u T ) 2 + (uf) 2 + |u| 2 a T + K| 2 uu r 
Hence 

Kit) = \J (\(£, u) 2 + i/i «£, u) 2 + |u| 2 |f| 2 )) ^ 

= i|p(a 2 (e,u) 2 + 6 2 (|u| 2 |e| 2 -(e,u) 2 ))^ = ^|p(« 2 (e,u) 2 + 6 2 |exu| 2 )^ 

R™ R™ 

= ^/p(« 2 (e^) 2 + & 2 Kxr / | 2 )^, 

R' 1 

where we denote 

. _. , ^ sin(a|£|£) A ., ,^, . „ sin(fo|£|t) „ 

^ = cos aft ¥>o + \'| 0i, 77 = cos 6 e * £o + ,\', £i, 

and for real vectors a = dei + a 2 e 2 + ... + a n e ni b = b i e 1 + 62^2 + ••• + b n e n , where 
{Gj}j = j-fi is the canonical basic of W 1 , we use the notation of the exterior product 



a x b = Y^ ( a i b j - aA)[ei x Gj ] G W 1 ^' 1 ^ 



Note that a x a = and we have the Lagrange's identity |a x b| 2 = |a| 2 |b| 2 — (a, b) 
Applying Plancherel theorem again for the kinetic energy functional, we get that 



U(t) = \ I r 



! a «> j ) 



L> 



d£. 
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Moreover, we have 

d _ tJ . s /_ sin(a|£|t) „ cos(a|£|t) „ \ ^ 



2 / sin (6|gt) cos(6|£|t) ,_ 

= («.<*>«+l«lV -(«,««, 
where 



?> ~ 6 iTi ^0 + TTTo Vi / ? 



lei 



iei J 



sin (a|£|t) „ cos(a|£|t) „ 
n = -a ^^ + ^P^, 



lei 



ki s 



sin^h , cos(6|e|t) ^ 



Note that (£, lei 2 /? - (£, /?> = 0. Hence 



iei s 



|-u(£,t) 

<9t v ^' ; 



= (e,«) 2 |e| 2 +||e| 2 /3-(e,/3)e| 
= (£, «) 2 |d 2 + ie| 4 |/?| 2 + & /3) 2 K| 2 - 2 <|e| 2 /?, (£, £> 

= iei 2 «e, «) 2 + iei 2 i/?r - & /?) 2 ) = lei 2 «e, «) 2 + I* x /?i 2 ) . 

Therefore, the total energy functional is calculated as 

E(t) = K(t) + C/(t) = \ fp(a 2 (£^> 2 + & 2 |£ x r# + |£| 2 «£, a) 2 + \£ x /3| 2 )) <£. 



Note that 

o 2 (e,^) 2 + le| 2 (e,«) 2 = (e,^) 2 + (e,|el«) 2 

/ <+> x ~ sin (a|£|t) „ 
£,acos(a|£|t)£ + ^^ 



ICI 



£, -asin(a|£|t)£ 



cos(a|£ |t) 

lei 



-V»] 



i=l 



E& ( acos(a|£|t)<£ ,i 



sin(a|£ |t) 

lei 

cos_(a]£jt) 

lei 



Vl, 



+ ( E & ( ~ a sin ( a l^l*) &M + 

= E ^ 2 (° 2 * + |Jp^?.i) + 2 E fcfc ( 



i ( a <fo,i<fo,j + 77j^^i,i^ij 



1 

IIP 



E ^^ 



and 

fc 2 |e x 77| 2 + |e| 2 |C x ^| 2 = le x 6r7| 2 H- |£x|0#| 

2 



.,_, . „ sin(fe|£|t) „ 
4>< (6cos(6|e|t)y + jj 1 ) <f 1 



i , . /,,^i \ - cos (6|£|t) „ 
0< ( -6sin(6|e|0v>o + ^p^i 
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+E (& (-&sin(6|e|t)^ 0j + C ° S £ ^.i) " (-*>Bw{b\Z\t)0 o ,i + 



01,i i 



cos(6|£|t) 

... , . ~^r~ 

So we obtain that 

E(t) = lJ P {a 2 r£^ A +^n2ti<piA + 

+ E f 62 te^o,j - £j<£cm) 2 + j7jl(&£i,i - £i&.i) 2 J ) ^ 

= 2 / p fl2 ( E &&w ) + E b2 (^oj - £m,if + E $.* ) d ^ 

ran \ \i=l / i<j 1=1 / 



2 

1 



P I a 2 



i=l % / «<i J 7 i=l / 



a 



Corollary 2.3. Assume that cp , cp 1: and f are vectors of Schwartz functions. Then, the so- 
lution to the elastodynamics wave equation (II .4P w?#/i fe zm'fz'fl/ condition (II .5P is unique. 

Proof. Suppose that u and u' are two solutions of (11.41 - I 1.5j) . It is clear that w = 
u — u' is a solution of the following homogeneous equation 

|^ - b 2 A\ w - (a 2 - 6 2 )Vdiv(w) = 0, 



w M)l i= o = -qT^^) 



= 0. 
t=o 



Applying Theorem 2.2., we obtain that the energy of w is equal to E(0) = 0. It 
implies that w = a.e., so we have completed the proof. □ 
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